errands$number.of.persons.by.home <- as.factor(errands$number.of.persons.by.home)
colnames(errands)[11]="number.of.errands.by.month"
colnames(errands)[14]="possess.vegetable.garden"
colnames(errands)[15:16]=c("participate.in.recycling.rainwater", "participate.in.waste.recycling")
colnames(errands)[29:34]=c("recognize.label.AB", "recognize.label.AOC", "recognize.label.AOP", "recognize.label.IGP", "recognize.label.LR", "recognize.no.label")
colnames(errands)[35:45]=c("buying.labelled.cheese", "buying.labelled.white.meat", "buying.labelled.red.meat", "buying.labelled.fish", "buying.labelled.alcoholic.drink", "buying.labelled.soft.drink", "buying.labelled.fruits.and.vegetables", "buying.labelled.chocolate", "buying.labelled.eggs", "buying.labelled.dairy.products", "buying.labelled.coffee.tea")

library(FactoMineR)
library(dynGraph)
library(EnQuireR)

summary(errands)

#MCA
errands.MCA<-errands[, c("buy.cheese", "buy.white.meat", "buy.red.meat", 
  "buy.fish", "buy.alcoholic.drinks", "buy.soft.drinks", "buy.chocolate", 
  "buy.fruit.and.vegetables", "buy.eggs", "buy.dairy.products", 
  "buy.coffee.tea", "label.decisive.criterion", "recognize.label.AB", 
  "recognize.label.AOC", "recognize.label.AOP", "recognize.label.IGP", 
  "recognize.label.LR", "recognize.no.label", "buying.labelled.cheese", 
  "buying.labelled.white.meat", "buying.labelled.red.meat", 
  "buying.labelled.fish", "buying.labelled.alcoholic.drink", 
  "buying.labelled.soft.drink", "buying.labelled.fruits.and.vegetables", 
  "buying.labelled.chocolate", "buying.labelled.eggs", 
  "buying.labelled.dairy.products", "buying.labelled.coffee.tea", 
  "no.labelled.product", "label.a.quality.product", 
  "label.certainty.on.the.origin", "label.a.specification", 
  "label.a.lifestyle", "label.a.preservation.of.the.gastronomy", 
  "label.local.products", "label.a.tradition", "label.a.better.taste", 
  "price.guarantee.of.quality", "price.margin.for.the.industrialist", 
  "price.earnings.for.the.producer", "price.emphasis.of.the.label", 
  "price.unjustified.superiority", "labelled.products.as.part.of.daily.diet", 
  "labelled.products.as.part.of.reception.of.guests", 
  "labelled.products.as.part.of.promotional.offers", 
  "labelled.products.as.part.of.discovery", "not.buying.labelled.products", 
  "sex", "age.bracket", "profession", "number.of.persons.by.home", 
  "number.of.inhabitants.in.the.city", "errands.in.hypermarket", 
  "errands.in.hard.discount", "errands.in.supermarket", 
  "errands.in.specialized.stores", "errands.in.market", 
  "number.of.errands.by.month", "limited.budget", "amount.of.the.errands", 
  "possess.vegetable.garden", "participate.in.recycling.rainwater", 
  "participate.in.waste.recycling")]

res<-MCA(errands.MCA, ncp=5, quali.sup=49:64, graph = FALSE)
plot.MCA(res, axes=c(1, 2), col.ind="#0000ff", col.ind.sup="blue", 
  col.var="#ff0000", col.quali.sup="#ff00ff", label=c("ind", "ind.sup", 
  "quali.sup", "var"), invisible=c(""), title="")
plot.MCA(res, axes=c(1, 2), choix="var", col.var="#ff8000", 
  col.quali.sup="#ce00ce", label=c("quali.sup", "var"), invisible=c(""), 
  title="")

plot.MCA(res, axes=c(3, 4), col.ind="#0000ff", col.ind.sup="blue", 
  col.var="#ff0000", col.quali.sup="#ff00ff", label=c("ind", "ind.sup", 
  "quali.sup", "var"), invisible=c(""), title="")
plot.MCA(res, axes=c(3, 4), choix="var", col.var="#ff8000", 
  col.quali.sup="#ce00ce", label=c("quali.sup", "var"), invisible=c(""), 
  title="")

#Hierarchical clustering
res.hcpc<-HCPC(res)

#dynGraph
dynGraph(res)

#EnQuireR
ENMCA(errands.MCA, quali.sup=50:64, report=TRUE)